Optical caustics of Kerr spacetime: the full structure 
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We present an exhaustive numerical investigation of the optical caustics in gravitational lensing 
by a spinning black hole for an observer at infinity. Besides the primary caustic, we examine higher 
order caustics, formed by photons performing one or several loops around the black hole. Our 
investigation covers the whole parameter space, including the black hole spin, its inclination with 
respect to the line of sight, the source distance, and the caustic order. By comparing our results with 
the available analytical approximations, we find perfect agreement in their respective domains of 
validity. We then prove that all caustics maintain their shape (a tube with astroidal cross-section) in 
the entire parameter space without suffering any transitions to different caustic shapes. For nearly 
extremal spin, however, higher order caustics grow so large that their cross-sections at fixed radii 
wind several times around the black hole. As a consequence, for each caustic order, the number of 
images ranges from 2 to 2(n + 1), where n is the number of loops spanned by the caustic. As for 
the critical curves, we note that for high values of the spin they develop a small dip on the side 
corresponding to prograde orbits. 

PACS numbers: 95.30.Sf, 04.70.Bw, 98.62.Sb 
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I. INTRODUCTION 



If General Relativity is the correct theory of gravity, 
the Kerr solution describes the spacetime metric outside 
spinning black holes [l[ . Therefore it is currently utilized 
in all models trying to reconstruct the phenomena ob- 
served around observed astrophysical black holes, either 
remnants of stellar collapse or supermassive black holes 
lying in the central regions of several galaxies 0, Q • A 
crucial step in the comprehension of physics in such ex- 
treme environments is the complete understanding of the 
phenomenology related to the bending of photon trajec- 
tories caused by spacetime curvature. In order to get 
reliable predictions on any observables, it is necessary to 
keep in mind that a Kerr black hole acts as a very strong 
gravitational lens, generating an infinite number of im- 
ages of the same source [||. The total flux of a source 
close to a Kerr black hole, (such as the accretion disk 
itself [5, 6] , an isolated bright spot on it @, H] or a star 
orbiting the black hole @ ) gets a significant contribution 
from the secondary image and from higher order images 
Q. This is also true for the details of fine structures 
in the profile of spectral lines, such as the iron Ka line 
in the X-ray domain, which are strong indicators of the 
presence of an intrinsic angular momentum of the black 
hole [hJ. 

On the other hand, the progresses in radio [ll| 
and infrared band [l2j interferometry, along with the 
projects of X-ray interferometry in space (MAXIM, 
http://maxim.gsfc.nasa.gov), foreshadow that resolved 
pictures of the nearest supermassive black hole (Sgr A* 
in the center of the Milky Way) will be feasible in a not so 
far future [l3[ . This will represent a spectacular advance 
in our knowledge of black hole physics. In particular, 
the contributions of different images of the same source 



will be identified and studied separately. Higher order 
images will then provide a huge amount of independent 
information on the inner portions of the accretion disks of 
supermassive black holes and will become precious wit- 
nesses of the strong gravitational field just outside the 
horizon. 

Gravitational lensing theory states that the multiplic- 
ity of the images of a given source depends on the source- 
lens-observer configuration. In a given spacetime met- 
ric, for an observer in a particular spacetime point, the 
multiplicity only depends on the source position. If the 
metric is stationary (as in the case of the Kerr metric) 
and the observer is static, the caustic can be defined 
in the 3-dimensional subspace at constant time as a 2- 
dimensional surface separating regions of space in which 
a source would give rise to a different number of images. 
When a point-like source crosses a caustic, a pair of ad- 
ditional images with infinite magnification is created or 
destroyed (the finite source size acts as a cut-off for real 
sources) |14| . 

It can be easily guessed that the study of the shape of 
the caustics is of fundamental importance for a reliable 
and complete description of the whole phenomenology 
related to the environment of astrophysical black holes. 
In fact, the multiplicity of the images and their brightness 
is essentially determined by the position of the source 
within the caustic structure. Even temporal variations 
in the observed overall luminosity may be due to caustic 
crossing of bright features around the black hole [l5[ • 

Surprisingly, 45 years after the discovery of the Kerr 
metric, the complete structure of the caustics of a spin- 
ning black hole has not yet been derived. Indeed, the 
complexity of the metric prevents from finding simple 
analytical solutions for the caustic surfaces. Even nu- 
merical studies are very challenging and not straight- 
forward. The first indication of the existence of non- 
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degenerate caustics came from the work of Cunningham 
and Bardeen who traced the light curves of a source 
star orbiting a spinning black hole. They noticed that 
the magnification of the primary and secondary images 
diverged at some particular points, signaling that caustic 
crossings were occurring. In spite of the huge number 
of ray-tracing codes in Kerr spacetime developed in so 
many years, the only comprehensive study of the caus- 
tic surfaces has been performed by Rauch and Bland- 
ford [TEl ]. They have explicitly shown that the primary 
caustic is a tube with a cross-section having the shape 
of an astroid (a closed curve with four cusps), which is 
very typical in gravitational lensing theory as soon as the 
spherical symmetry of the lens is broken by an external or 
internal perturbation. They have shown several pictures 
of the primary caustic and derived some simple asymp- 
totic behaviors for its size. Besides the primary caustic, 
the authors have mentioned the existence of higher order 
caustics, but they have not shown any pictures of them, 
leaving several questions about the size and the shape of 
these caustics open. Later on, Sereno and De Luca have 
found an analytical approximation for the primary caus- 
tic valid for large source distances [16[ . In a series of pa- 
pers based on the strong deflection limit approximation 
[17| . we have derived a perturbative analytical approxi- 
mation describing the higher order caustics (but not the 
primary caustic) [18l.ll9l.l20l]. These approximations show 
that at low spin values the higher order caustics are still 
tubes with astroidal cross sections with increasing size. 
Basically, this is all we know about caustics in the Kerr 
spacetime. 

All these studies provide several hints about the caus- 
tic structure of the Kerr black hole lens in several limits. 
Yet, the fate of the higher order caustics at high values of 
the spin still remains unclear. As they become larger and 
larger with increasing spin, do they undergo any transi- 
tions to different caustic shapes? Does their size stay 
finite? Do they merge? The present work provides clear 
answers to these and other questions of theoretical and 
observational relevance, clarifying the whole panorama of 
the caustic structure of the Kerr spacetime. We present 
a thorough numerical analysis of the caustics generated 
by a Kerr black hole at all caustic orders, studying their 
dependence on the black hole spin, its inclination and 
the source distance. The reliability of our results is also 
double- checked against former studies and all analytical 
approximations available up to now. 

The paper is organized as follows. Section II traces the 
methodology followed for the generation of the caustics, 
referring to the appendix for a detailed explanation of all 
steps. Section III deals with the primary caustic. Section 
IV discusses the dependence of higher order caustics on 
the spin and its inclination. Section V focuses on the 
caustics of extremal black holes. Section VI is devoted to 
critical curves in the observer's sky. Section VII contains 
the conclusions. 



II. METHODOLOGY 



The Kerr metric in Boyer-Lindquist coordinates [21[ is 
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M is the mass of the black hole and a is the specific 
angular momentum, ranging from (Schwarzschild black 
hole) to M (extremal black hole). The horizon radius is 
r h = M + ^M 2 - a 2 . 

Let us consider a static observer at coordinates 
(r Q , ■do, 4> ) with (j) a = and let us assume r a ^> 2M for 
simplicity. In practice, all astrophysical situations sat- 
isfy this limit. The observer's polar coordinate also 
coincides with the inclination of the spin axis with re- 
spect to the line of sight between the observer and the 
black hole. It is often convenient to work with the coor- 
dinate /i = cost? instead of 1). Therefore, one can define 

fig = COS^Q. 

Similarly, we consider a source at coordinates 
(r 8 ,i?«,0 s )) with Ps = cos?9 s . 

The observer can construct his own coordinates (0\, 62) 
in the sky, such that the black hole is in the origin and 
the 62 axis coincides with the projection of the spin of 
the black hole on the sky. 

Now we make a qualitative summary of the numerical 
algorithm for finding the caustics and state a few ba- 
sic equations that should be sufficient to understand the 
main points. In Appendix A we report all the details of 
the calculation and the definitions of all variables. 

The first important concept to be reviewed is the 
shadow of the black hole [2| . Imagine tracing a photon 
trajectory from the observer back toward the black hole. 
The trajectory is completely determined by the angles 9\ 
and 82, which specify the final direction of the photon 
when it hits the observer. These angles may be consid- 
ered as the initial conditions for tracing back the photon 
geodesic in the Kerr metric. Then, it is well-known that 
there exists a set of photon trajectories that are traced 
back to grazing unstable fixed radii orbits around the 
black hole. These particular trajectories define a locus in 
(01,82), which is called the shadow border of the black 
hole. As evident from Fig. [TJ in the limit a — > 0, this 
locus has a circular shape with angular radius 3\^3M/r . 
For non-vanishing spin, it is slightly displaced toward 
the right direction (remember that we have chosen coor- 
dinates such that the spin of the black hole is projected 
on the positive 8 2 axis) and becomes slightly flattened 
on the left side. The shape of the shadow depends on 
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FIG. 1: The shadow border as seen by an equatorial observer 
for different values of the black hole spin. 9i and 62 are the 
angular coordinates in the observer's sky in units of 2M/r . 
The black hole is in the origin and the spin axis points in the 
positive 62 direction. 



the black hole spin and its inclination on the line of sight 
and can be used, in principle, to extract these parameters 
@; EH H3| ■ The reason for the name "shadow border" is 
that all sources lying outside the unstable circular orbits 
generate images appearing outside this locus in the ob- 
server's sky. Only sources very close to the black hole 
may have images inside the shadow border. Therefore, 
to an instrument with the sufficient angular resolution, 
the black hole should appear as a dark spot (shadow) sur- 
rounded by a more luminous region [l3T ] . In the appendix, 
we review the derivation of the shape of the shadow bor- 
der, which is explicitly given by Eqs. (|A11[) and (|A12[) 
in parametric form as functions of r m , which represents 
the radius of the unstable photon orbit from which the 
photon emerges and hits the observer at angles (9\, 6*2). 

It is convenient to switch from r m to a more friendly 
parameter that we call 77 as defined by Eq. (|A14[) . 77 is 
very similar to an angular variable: as it ranges in [— ir, n], 
the whole shadow border is spanned. The intersections 
with the Q\ axis are at 77 — (with 9\ > 0) and 77 = ±ir 
(with 6*i < 0). rj = it/2 corresponds to a point close to 
the positive 62 axis and 77 = —n/2 corresponds to a point 
close to the negative 62 axis. As the spin axis is projected 
to the positive #2 axis, we can deduce that photons reach- 
ing the observer from the left of the black hole (77 ~ ±7r) 
come from orbits co-rotating with the black hole, whereas 
photons reaching the observer from the right of the black 
hole (77 ~ 0) come from orbits rotating in the opposite 
sense with respect to the black hole. We shall often re- 



fer to photons on co-rotating orbits as prograde photons 
and to photons on counter-rotating orbits as retrograde 
photons. Photons reaching the observer from above or 
below the black hole (77 ~ ±7r/2) come from quasi-polar 
orbits. 

The shadow border can be used to construct new co- 
ordinates (tp, rj) in the observer's sky, such that the sky 
coordinates {B\, 62) can be re-expressed in terms of them 



h = 0i(il>,rj) 
h = 2 (tp,v)- 



(4) 
(5) 



These coordinates are similar to polar coordinates, with 
77 playing t he role o f the polar angle and tp related to the 
modulus y/0f + 6*2 ■ The choice of the precise form of the 
radial variable tp as a function of \JQ\ + Q\ is done with 
the purpose of simplifying the numerical calculations, as 
will be explained in the following paragraphs. 

A detailed study of the geodesies equations at fixed 
source distance r s (see Appendix) provides the sought 
after relation between the source's angular coordinates 
and the variables ip an d 77 



<p s = <p s (lp,V)- 



(6) 
(7) 



As ip and 77 determine the position of the images in the 
observer's sky, Eqs. © and |(7J) represent the Kerr lens 
mapping. 

Of course, the detailed expression of this mapping may 
become very involved and cumbersome to numerical cal- 
culations if one does not make a proper choice of the vari- 
able ip. Luckily, the explicit form of the function (i s in 
Eq. ([5]) consists of an oscillating Jacobi elliptic function 
whose argument monotonically increases when yflf+5| 
decreases. From the physical point of view, this means 
that the closer the photon passes to the black hole the 
larger is the number of oscillations in the polar motion 
it performs. Then, it is natural to define tp as the argu- 
ment of this Jacobi elliptic function. With this definition, 
the function /j, s (tp,r)) resembles a sinusoid in tp, with the 
amplitude of the oscillations determined by 77. 

With a suitable normalization of tp, we can set the 
period of the Jacobi elliptic function to Aip = 2. In this 
way, we can establish the equivalence ji s = ip = m, 
with m integer. Furthermore, the number of inversions 
in the polar motion of the photon from the source to the 
observer is m = [tp + 1/2], where by [x] we indicate the 
integer part of x. The precise definition of tp is given 
by Eq. (|A36[) . though it is necessary to go through a 
conspicuous portion of the appendix to understand all 
the details. 

As long as the Jacobian determinant of the lens map- 
ping is different from zero, the mapping is locally invert- 
ible and the number of images stays constant. Images 
can be created or destroyed only in the critical points 
where the Jacobian determinant vanishes. It turns out 
that for each value of m there is one Jacobian critical 
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FIG. 2: Cross sections of the primary caustic for a — 0.9998M 
and \i = 0. From left to right, the radial coordinate is r 3 = 
5M, 4.5M, AM, 3M, 2.3M, 1.82M, 1.42M, 1.22M. 



point in the interval between two consecutive polar in- 
versions to — 1/2 < if) < m + 1/2, for all values of r\. 
Therefore, we can identify to with the caustic order. The 
primary caustic is obtained with a single polar inversion 
(to = 1), the second order caustic is generated by photon 
trajectories with two polar inversions (to = 2), and so on. 
With the critical points found in the (ip, rf) space, we can 
readily draw the critical curves in the observer's sky by 
applying Eqs. (H]) and ([5]), and the caustics by applying 
Eqs. © and 0. 

In conclusion, once we fix all the parameters (black 
hole spin a, spin inclination /!„, source distance r s , caus- 
tic order m), we are able to trace the corresponding caus- 
tic and critical curve through suitable numerical calcula- 
tions. 



III. PRIMARY CAUSTIC 

The primary caustic is the caustic generated by light 
rays with a single inversion point in the polar motion. 
In the Schwarzschild limit (a — > 0), it degenerates to a 
line starting from the black hole and extending to infin- 
ity in the opposite direction with respect to the observer 
(cf) = —if). Distant sources that approach the primary 
caustic experience the standard weak deflection gravita- 
tional lensing as described in classical textbooks. The 
two weak deflection images merge when the source lies 
right on the primary caustic, forming the well-known Ein- 
stein ring. 

If we turn the black hole spin on, the degenerate one- 
dimensional caustic becomes a finite-thickness tube with 
a cross-section having the shape of a four-cusped astroid. 
Several cross-sections at fixed radii are shown in Fig. 
[5] for a nearly extremal Kerr black hole. If the source 
lies inside the enclosed region, two additional images are 
present 

As already shown in the extensive study by Rauch & 
Blandford [15j . the cross-section is reflection-symmetric 



FIG. 3: The primary caustic tube for a = 0.9998M and [i = 
0. In building this 3-dimensional representation, we have cho- 
sen coordinates such that x = r sin i? cos <j>, V = rsini9sin0, 
z = roostf, following Ref. [l5| . The spin axis is directed to- 
ward the top. The straight solid line indicates the direction 
towards the observer, whereas the dashed line points in the 
opposite direction. The primary caustic for a Schwarzschild 
black hole (a — 0) coincides with this dashed line. 



at very large distances, but becomes more and more dis- 
torted approaching the black hole. Moreover, the caustic 
tube is shifted clockwise (i.e. in the sense opposite to the 
rotation of the black hole), as evident from Fig. [2 where 
it can be noted that the caustic does not lie at <j> = —it, as 
in the Schwarzschild case (We recall that we have put the 
observer in <f> = 0). This shift becomes larger and larger 
for sources closer to the horizon. A full 3-dimcnsional 
picture of the primary caustic tube was first presented 
in Ref. [l5| and is shown again here in Fig. [3J It can 
be noted that the caustic tube becomes belt-like as it 
winds around the black hole horizon, represented by the 
spherical surface in the figure. 

As anticipated in the introduction, the caustic sepa- 
rates two regions of space in which a source gives rise to 
a different number of images. After drawing the primary 
caustic, it is useful to explain the way it affects the for- 
mation of the images around the corresponding critical 
curve. We remind the reader that only the images of 
order to = 1 (i.e. generated by photons with one polar 
inversion) are involved. 

In Fig. [4] we have drawn a schematic picture of the 
images of a source close to the caustic as detected by 
an observer. In Fig. 0^ we have a cross-section of the 
primary caustic at r s = 12M, along with several possible 
source positions, represented by different symbols. In 
Fig. UJd we show the corresponding critical curve as a 
solid line and the shadow border in dashed style. The 
images corresponding to the source positions of Fig. 0^ 
are shown with the same symbol. We also remind the 
reader that the spin is directed toward the north pole 
(/j = 1) in Fig. QJi and is parallel to the 62 axis in Fig. 
Sb. 

If the source is outside the caustic (the cross in Fig. QJl) 
there are two images (the two crosses in Fig. HJd). One 
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FIG. 4: A schematic picture of the images of a point-source 
at r s — 12M generated by a black hole with spin a = 0.2M as 
seen by an equatorial observer, (a) The primary caustic cross- 
section and four possible positions of the source, indicated by 
the cross, the filled circle, the empty circle and the box. (b) 
The corresponding images in the observer's sky together with 
the critical curve (solid line) and the shadow of the black hole 
(dashed line). 



image is on the same side of the source and lies outside 
the critical curve (primary image) and the other is oppo- 
site and lies inside the critical curve (secondary image) . If 
the source enters the caustic from the right (empty box), 
two images with opposite parities are created on the left 
of the critical curve (there is a reflection with respect to a 
vertical axis from the caustic to the critical curve). If we 
place a source on the equatorial plane close to the right 
cusp (filled circle), we have three images close to the left 
tip of the critical curve, corresponding to prograde pho- 
tons. Finally, if we have a source close to the upper cusp 
(empty circle), we have three images close to the upper 
tip of the critical curve. 

Then, to summarize, the right cusp of the caustic in- 
volves the creation of images generated by prograde pho- 
tons; the opposite occurs for the left cusp. Therefore, 
we will often refer to the right cusp as the prograde cusp 
and the left cusp as the retrograde cusp. As for the upper 
and lower cusps, they are associated with the formation 
of images generated by quasi-polar orbits. In the pri- 
mary caustic (and in all odd higher order caustics), the 
upper cusp is associated with the upper tip of the critical 
curve, whereas in all even order caustics it is associated 
with the lower tip (there is also a reflection with respect 
to a horizontal axis in even order caustics [Hj]). 

This picture is a useful reference to understand how the 
different parts of the caustic are associated with the the 
corresponding parts of the critical curve and where the 
images are created or destroyed when the source crosses 
a caustic. From now on, we will concentrate mostly on 
the shape of the caustics, confining the discussion about 
the critical curves to Section IVT1 

Coming back to the primary caustic surface, we find 
that for any values of the black hole spin the primary 
caustic performs an infinite number of turns around the 
horizon. This fact is evident from Fig. [5^i, where the 
azimuthal coordinate of the left cusp (corresponding to 
retrograde photons) is plotted vs the logarithm of the 
difference between the radial coordinate and the horizon 
radius r/j. This plot shows that at large distances the 
primary caustic always tends to — 7r, whatever the value 
of the spin. In the opposite limit, very close to the black 
hole horizon, the position of the retrograde cusp shifts 
clockwise more and more without stopping ((f) increases 
with a negative sign). We can easily establish a loga- 
rithmic law 4> c ~ c(a) log(r s — r/j), where the coefficient 
c(a) depends on the black hole spin, becoming larger and 
larger for high values of the spin. 

As already guessed in Ref . [26] , the divergence of the 
azimuthal position of the caustic in the approach to the 
horizon is a direct consequence of the divergence of the 
radial integral I2 (Eq. (|A33[) in the appendix) appearing 
in the lens equation for <fi s (|A38[) . This integral contains 
a factor 1/A, which generates the logarithmic divergence 
of I2 ■ This has been explicitly shown here by numerical 
calculations for the first time. 

In Fig. [SJd we plot the size of the primary caustic 
(calculated as the distance between the right and the left 
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FIG. 5: (a) Azimuthal coordinate of the retrograde cusp of 
the primary caustic as a function of the radial coordinate, 
(b) Size of the primary caustic as a function of the radial 
coordinate. In both figures, from bottom to top, the curves 
are drawn for a = 0.2M, 0AM, 0.6M, 0.8M, 0.9998M. 
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FIG. 6: Comparison between numerical calculations (solid 
lines) and the analytical approximation by Sereno and De 
Luca [jij (dashed lines), (a) Position of the retrograde cusp 
as a function of the radial coordinate, (b) Size of the caustic 
as a function of the radial coordinate. From bottom to top, 
the curves are for a = 0.2M, 0AM, 0.6M, 0.8M, 0.9998M. 



cusp) as a function of the radial coordinate in logarithmic 
scale. We see that the size tends to zero at very large 
distances, whatever the value of the spin. This is only 
true for the primary caustic, whereas all higher order 
caustics tend to a finite size at large distances, as we 
shall see in the next section. Close to the horizon, the 
size of the primary caustic tends to a constant. 

Recently, an interesting analytical approximation for 
the primary caustic has been proposed by Sereno and 
De Luca in Ref. [la ] . using methods previously devel- 
oped in Ref. [23| . This approximation is obtained by an 
expansion in the parameter 



caustic reads 



£ = 



>M(r 



(8) 



which is valid at large distances from the black hole for 
any values of the spin. The explicit expression of the 
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with ?7 ranging from — 7r to n. 

In Fig. [5] we compare our numerical results with the 
formulae by Sereno and De Luca, finding a perfect agree- 
ment at large distances from the black hole. In partic- 
ular, Fig. [6^ represents the azimuthal position of the 
retrograde cusp as a function of the radial coordinate in 
a log-log plot, whereas Fig. [6Jd shows the size of the caus- 
tic as a function of the radial coordinate in a log- log plot. 
The approximation only fails at very small distances from 
the horizon. It is interesting to note that the exact size 
is smaller than that predicted by Eq. © for small spins, 
whereas the situation is reversed at higher spins. 
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IV. HIGHER ORDER CAUSTICS 

It is well-known that a Schwarzschild black hole gen- 
erates two infinite sequences of lensed images of a given 
source [24j . Apart from the two main images, usually 
treated in the Weak Deflection Limit, all remaining im- 
ages can be described in the Strong Deflection Limit 
@> G3, H3] • Higher order images arise from the fact that 
photons grazing the unstable photon orbit may perform 
one or more turns around the black hole before emerg- 
ing. For each turn, we have an additional pair of images, 
whose luminosity decreases exponentially as we increase 
the number of turns. In practice, gravitational lensing 
phenomenology is replicated for each number of turns 
around the black hole, including a new ring-like criti- 
cal curve and a degenerate caustic line. Furthermore, 
the structure is replicated even for sources in front of 
the black hole (this geometrical configuration has been 
dubbed retro- lensing [25j]). In this case, apart from the 
direct image of the source, all other images are due to 
strongly deflected photons turning around the black hole 
one or more times and reaching back the observer. 

Standard lensing and retro-lensing caustics can be 
treated in a unified way, noting that we have a new 
caustic each time we add a new inversion point in the 
polar motion of the photon. As anticipated in Section 
II, the number of inversions can then be identified with 
the caustic order. For example, the primary caustic is 
generated by photons with a single inversion point, the 
second order caustic (first retro-lensing caustic) is gen- 
erated by photons with two inversion points, the third 
order caustic (second standard lensing caustic) is gener- 
ated by photons with three inversion points, and so on. 
This situation holds even when we switch the angular 
momentum on. Nevertheless, analogously to what hap- 
pens for the primary caustic, higher order caustic tubes 
are shifted and acquire a finite thickness. Up to now, 
higher order caustics have been studied in the Strong De- 
flection Limit approximation [26] in the equatorial plane 
and successively for small values of the black hole spin 
in a fully analytical way [TH, [l^, [2(| • In this section, we 
shall present a complete study of higher order caustics in 
the whole parameter space. 

First we note that at large distances there is a funda- 
mental difference between the primary caustic and the 
higher order caustics. In fact, whereas the primary caus- 
tic shrinks to zero size and tends to the Schwarzschild 
point-caustic at (f> = — n for r s — > oo, each higher or- 
der caustic tends to a well-determined asymptotic astroid 
shape in the coordinates (</j, [i). Therefore, rather than a 
tube, at very large distances the 3-dimensional region en- 
closed by the caustic resembles a pyramid with astroidal 
base and vertex in the black hole. The volume enclosed 
within the caustic grows as O c r^ , where tt c is the angular 
area within the asymptotic astroidal shape. 

Let us begin by studying the dependence of the asymp- 
totic cross-sections of the caustics at very large distances 
on the black hole spin a. In Fig. [7] we show the second 
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FIG. 7: Second order (a) and third order asymptotic 
caustics (b) for an equatorial observer fi — and dif- 
ferent values of the spin. From left to right, a — 
0.02M, 0.2M, AM, 0.6M, 0.8M, 0.9998M. 



order caustic (first retro-lensing caustic) and the third 
order caustic (second standard lensing caustic) for differ- 
ent values of the spin. At very low values of a, higher 
order caustics are very small and only slightly displaced 
from the Schwarzschild positions 4> = —mir. As a grows, 
the size and the displacement of the caustics grow. At 
the same time, the caustics become asymmetric, with 
the cusp on the right (corresponding to prograde pho- 
tons) more stretched than the cusp on the left (retrograde 
photons) . 

At nearly extremal spins, higher order caustics become 
very large but do not undergo any transitions to different 
shapes. This regime is particularly interesting and will 
be deeply investigated in the following section. 

Now we turn to the dependence of the asymptotic caus- 
tics on the spin inclination with respect to the line of 
sight. In Fig. [5] we show some asymptotic cross-sections 
of the caustics for different values of the spin inclina- 
tion \i Q . The retro-lensing caustics are always centered 
on i.e., they lie at the same latitude as the observer, 
whereas the standard lensing caustics lie at opposite lat- 
itude. Both retro-lensing and standard lensing caustics 
preserve the same shift in the azimuthal direction as the 
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FIG. 8: Second order (a) and third order asymptotic caustics 
(b) generated by a Kerr black hole with spin a = 0.6M for 
different values of the inclination fi = cos # . From the center 
to the periphery, fi a = 0, 0.3, 0.6, 0.9, 0.99. 



Ho = caustic. Increasing /j from (observer on the 
equatorial plane) to 1 (observer on the polar axis), all 
higher order caustics gradually shrink to zero size. An ob- 
server on the polar axis would thus have point-like caus- 
tics on the optical axis as in the Schwarzschild case. An- 
other interesting feature is that the caustics are slightly 
stretched on the side closer to the equatorial plane. This 
is not just an effect of the use of the projected coordinate 
fi as it can be seen even in terms of the polar angle 
We must conclude that multiple images arise more easily 
for sources close to the equator. 

In order to double-check our numerical calculations, we 
can compare our asymptotic caustics with some approx- 
imation schemes. The first is the one developed in Ref. 
[26| . based on the Strong Deflection Limit. In that paper, 
only the intersections of the caustics with the equatorial 
plane have been calculated numerically. Nevertheless, a 
comparison with these results allows us to check the ret- 
rograde and prograde cusp positions with a completely 
independent calculation. The second approximation is 
provided by the fully analytical formulae derived in Refs. 
13 [H, H3] in the Strong Deflection Limit with an ex- 
pansion to second order in the black hole spin a. They 



H. = (-l) m L + Rm(l - Mo 2 ) 3/2 sin 3 r, 



(11) 



4> s = -mir - Acj) m + Rms/l - /i Q 2 cos 3 77, (12) 
where the azimuthal shift is 

A0 m =_|^| + 21og(2V3-3 

. \ (2y/r; + y/6M + r s 
log — - — = — 
6 [ 3^/(r s - 2M) 

and the semi-amplitude of the caustic is 

R rn = a 2 ^-(5tott + 8V3-36) 
18 

(9M + 2r s - 



(13) 



3\/VrW6M + r s 



(14) 



In the limit r s ^S> 2M, these expressions reduce to their 
first lines respectively. 

Fig. [5] shows the dependence of the positions of the 
cusps of the asymptotic caustic on the black hole spin. 
Along with our numerical results, we have plotted the 
two approximations described above. We see that the 
Schwarzschild limit is correctly reproduced as a — > 0, 
since the caustics become point-like and return to the 
optical axis (— 2-7T for the second order caustic and — 37r 
for the third order caustic). As soon as we switch the 
spin on, the caustics are shifted clockwise (more negative 
4>) and acquire finite extension, since the prograde cusp 
is more shifted than the retrograde cusp. 

The perturbative approximation works quite well up to 
a ~ 0.2M. At higher spins, it underestimates the shift of 
the caustics. On the other hand, the SDL approximation 
without the expansion in powers of a works extremely 
well for the third order caustic, as the solid and dotted 
lines are practically indistinguishable in Fig. [5Jd. For 
the second order caustic, instead, it works quite well for 
the retrograde cusp but fails at moderate and high spin 
values for the prograde cusp. We must keep in mind 
that the SDL approximation is designed for high caustic 
orders and is just marginally applicable to the second 
order caustic. 

Now let us come to the full development of the caustic 
surfaces, from very far source distances up to the hori- 
zon. In Fig. [TUlwe show the second order caustic tube 
for a rather small value of the spin (a = 0.2M). It looks 
quite similar to the primary caustic. The caustic tube 
has astroidal cross-section and winds an infinite number 
of times around the horizon. However, the cross-section 
does not shrink to zero size at infinity, but preserves the 
same angular size. A 3-dimensional picture containing 
several caustics at the same time would not be so read- 
able, as all caustics wind around the horizon so tightly 
that they would be undistinguishable. However, a picture 
giving the idea of how several caustics wind around the 
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FIG. 9: Azimuthal position of the prograde and retrograde 
cusps of the asymptotic caustic as a function of the black hole 
spin, (a) Second order caustic; (b) Third order caustic. The 
solid lines are from numerical results in this paper, the dotted 
lines are the SDL approximation of Ref. [2^ |. the dashed lines 
are the perturbative results of Ref. . 



FIG. 10: The second order caustic surface for a = 0.2M and 
/i = 0. The straight line indicates the direction towards the 
observer. 




black hole is in Fig. QT] where we have represented the 
projection of the caustics on the equatorial plane. The 
radial coordinate has been replaced by the logarithmic 
coordinate log(r s — rh)/2M + 11 in order to put the "spi- 
ral arms" of the caustics in better evidence. Of course, 
such a representation becomes meaningless at distances 
less than r/, + 10 _11 M. We have represented the pri- 
mary caustic (labelled by "1") along with higher order 
caustics up to the sixth order. In this way we have explic- 
itly shown where the caustics lie with respect to the ob- 
server's direction (labelled by "O"). In the Schwarzschild 
limit, all caustics reduce to a line: odd ones are on the 
opposite side and even ones (retro-lensing caustics) are 
on the same side as the observer. At a > 0, all caus- 
tics are shifted clockwise and acquire a finite extension, 
proportionally to their order. 

Let us come to a more quantitative analysis. Fig. 1121 
shows the azimuthal position of the retrograde cusp as 
a function of the radial coordinate for the second and 
third order caustic. The behavior is very similar to the 
primary caustic, with the logarithmic divergence for r s — > 
r/j. However, at large radii, the retrograde cusp does not 



FIG. 11: A representation of the caustics on the equatorial 
plane in which the radial coordinate has been replaced by 
log(r s — rh)/2M + 11. The outer circle represents the sur- 
face at log(r s — rh)/2M = 4. Inner circle are in steps of 1 
in such a logarithmic coordinate. The observer is in the di- 
rection labelled by O on the equatorial plane. The number 
at the end of each caustic represents the respective caustic 
order. Here the spin is a — 0.2M with the black hole rotating 
counterclockwise. 



settle to the Schwarzschild position — rn7r, but a finite 
shift remains. In Fig. [T3]we plot the size of the caustics 
(estimated as the distance between the prograde and the 
retrograde cusp) as a function of the radial coordinate. 
Indeed, the size of the caustics tends to two different 
constants at large and small radii. This behavior has 
already been noted for small a caustics in Ref. [20| . 

We conclude this section by a picture representing the 
relative error in the prediction of the position and size of 
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the caustics by the perturbative formulae (flU)) and (fH)) . 
In Fig. [14] we can see that the error decreases expo- 
nentially with the caustic order m, as predicted by the 
SDL approximation. However, a residual error remains 
because of the expansion in powers of a. In any case, 
for the value a = 0.02M used in this plot, the error is at 
most of 3% for the second order caustic at large distances. 
Therefore, the SDL approximation works amazingly well 
at all radial distances, reproducing the logarithmic spi- 
ralling perfectly. 



V. CAUSTICS OF EXTREMAL KERR BLACK 
HOLES 

As shown in the previous sections, higher order caus- 
tics become very large at high values of the spin. Since 
there are several indications that nearly extremal spin- 
ning black holes may be not uncommon in the universe, it 
is particularly important to focus on the caustic structure 
of this particular limiting species of Kerr black holes. 

First, we show the asymptotic cross-sections at large 



source distances for several caustic orders. Apart from 
the second order caustic, all higher order caustics extend 
for more than a complete loop around the black hole. 
The left (retrograde cusp) is shifted clockwise from the 
Schwarzschild position linearly with the caustic order to, 
following the law 

0/ = -3.86m + 0.87. (15) 

On the contrary, the right (prograde) cusp is shifted ex- 
ponentially with to, following the law 

<t> r = -1.05m - 0.36 exp(1.57m). (16) 

Caustics overlap each other without any interference, as 
they correspond to trajectories with different numbers of 
inversion points in the polar motion, which live in dif- 
ferent regions of the (ip,rj) space. Also the extension in 
latitude increases, with the northern and southern cusps 
gradually approaching the poles. 

The reader may wonder about the physical conse- 
quences of having caustics extending for more than 2tt 
in azimuth. In Fig. [TB] we plot the third order caustic 
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FIG. 15: Asymptotic caustics of an extremal Kerr black hole 
for an equatorial observer. From left to right the caustic order 
is m = 2, 3, 4, 5. 
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FIG. 16: Third order asymptotic caustic of an extremal Kerr 
black hole for an equatorial observer in the <j>-[i plane. In this 
plot, we have taken into account the fact that the azimuthal 
coordinate is periodic with period 2n. 

reporting all points in the physical interval < <\> < 2tt. 
In this way, the caustic appears to exit from the right 
side and re-enter from the left side for five times before it 
ends in the prograde cusp. Suppose to have a source on 
the same side of the observer (^> s ~ 0) and close enough 
to the northern pole (fj, s ~ 1) to be outside of the caus- 
tic surface. In this situation, there are only two images 
formed by photons with three inversions in the polar mo- 
tion. If we move the source down, by decreasing ^i s , the 
source meets the caustic for the first time. At this point, 
two new images appear, as we know from gravitational 
lensing theory. Their initial position on the critical curve 
is the one that corresponds to the caustic point where 
the source crosses the caustic. However, if we continue 
to move the source down, at some point we meet the 
caustic for the second time. Two new images appear in 
a different position on the critical curve. In particular, 
these new images are still formed by photons with three 
inversion points in the polar motion but they perform 
one more loop in the azimuthal motion before reaching 
the observer with respect to the previous pair. If we 
continue to decrease /i s , we can cross the caustic 4 more 
times, including the final cusp crossing. Therefore, for an 
extremal Kerr black hole, we can have up to 14 images 
with three inversion points in the polar motion. 

The number of times a particular caustic can be 
crossed is simply given by the number of azimuthal loops 
spanned by the caustic plus one. With the empirical for- 
mulae fll~5j) and (HHJ), we can estimate that a caustic of 
order m may generate up to 

n = [0.36 exp(1.57m) - 2.82m + 0.87] (17) 

new pairs of images. In practice, for a source very close 
to the pole, the number of images with m inversions in 
the polar motion will always be the minimal one, i.e. 2. 
On the other hand, for a source close to the equatorial 
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FIG. 17: Second order (a) and third order (b) asymptotic 
caustics of an extremal Kerr black hole for different positions 
of the observer. From the center to the periphery the caustics 
are drawn for fi = 0, 0.3, 0.6, 0.9, 0.99. 




plane, the number of images will always be very close 
to the maximal one 2(n + 1). It is interesting to connect 
this result with the exponential decrease in the magnifica- 
tion of higher order images with the caustic order. From 
one side we have an exponential increase of the number 
of images. From the other side we have an exponential 
decrease in the magnification factor of each image. How- 
ever, even if we have not calculated the details of the 
magnification, we can expect that the decrease in the 
magnification dominates on the increase in the number 
of images, in such a way that the total flux in the images 
of order m still decreases exponentially with the caustic 
order. Finally, we also note that, apart from the first 
pair of additional images, all images created by repeated 
caustic crossings occur in a region close to the stretched 
prograde cusp. Therefore, all these additional images will 
appear to the observer on the left side of the black hole. 

As for the dependence of the caustics on the spin incli- 
nation on the line of sight, we have plotted several second 
order and third order caustics in Fig. Q~7] As in Fig. [5J 




FIG. 18: 3-dimensional pictures of the second order caustic 
surface of an extremal Kerr black hole for an observer on 
the equatorial plane. The surface has been plotted for radial 
distances in the range [2.2M, 20M]. (a) View from the north 
pole, (b) View from the observer side, lifted by 10° above the 
equator, (c) View from the left side. The straight line in (a) 
and (c) points towards the observer. 



the caustics keep the same azimuthal shift for their cen- 
tral region while their size decreases with increasing fi a . 
Even if the fi a — caustics extend for several azimuthal 
loops, by increasing fi a up to fi a = 1, all caustics shrink 
to a point. Furthermore, we can see in particular from 
Fig. 117b that the caustics maintain their extension in lat- 
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FIG. 19: Azimuthal position of the retrograde cusp (a) and 
prograde cusp (b) as a function of the logarithm of the dis- 
tance from the event horizon. From bottom to top, the curves 
are for caustic order m = 2,3, 4, 5. 



itude until the azimuthal stretch is reduced to less than 
one loop. Only then, the latitudinal size starts to de- 
crease significantly. As a consequence, since the left and 
right cusps are always at fj, = ±/x OJ we have a consider- 
able asymmetry between the side of the caustic towards 
the pole and the side towards the equator, which remains 
much larger. 

After the study of the asymptotic cross-section of the 
caustic surfaces, it is time to examine the full develop- 
ment of the caustics of an extremal Kerr black hole. As 
it can be deduced from the extension of the asymptotic 
cross-sections, the caustic surfaces become so wide that 
they no longer resemble a tube. They assume an almost 
disk-like shape with the thickness increasing with the dis- 
tance from the black hole. Some 3-dimcnsional pictures 
of the second order caustic are shown in Fig. [TH] with 
different points of view. The caustic has been calculated 
for an observer on the equatorial plane (whose direction 
is indicated by the straight line). One can appreciate 
the half-disk shape which winds around the black hole 
more and more. The retrograde cusp is the equatorial 



cusp on the left side, which also enjoys the wider lat- 
itudinal extension. The prograde cusp is on the right 
side, where the caustic has a spiky cross-section. The 
prograde cusp winds around the black hole much more 
quickly than the retrograde cusp, at least initially. It can 
be easily imagined that the pictures of even higher order 
caustic surfaces, which spread over several loops around 
the black hole, are very difficult to show in a clear way. 
The caustic surface self-intersects several times. 

Finally, we plot the position of the retrograde and pro- 
grade cusps in Fig. Q1J] We can see that both retrograde 
and prograde cusps wind around the black hole following 
the usual logarithmic law (f> ~ clog(r s /r/j — 1). However, 
the prograde cusps start to follow the logarithmic law at 
much larger distances than the retrograde cusps. This 
means that the caustic has a very large increase from the 
asymptotic size at large radial distances to the final size 
at distances very close to the black hole. The slope of 
the logarithmic laws is the same for all prograde and ret- 
rograde cusps and amounts to c = 57.55, which means 
that each time we decrease the distance from the horizon 
by a factor 10, the caustic makes 9.16 turns around the 
black hole. 



VI. CRITICAL CURVES 

As stated in Section II, once we find a critical point of 
the lens mapping in the (ip, 77) space, we can easily de- 
termine the corresponding critical point in the observer's 
sky, spanned by the angular coordinates (#i,02)- As rj 
varies from — ir to 71", the critical points describe a closed 
curve around the black hole. Therefore, for each caustic 
order m, for a given source distance r s , we have a distinct 
critical curve in the observer's sky. In particular, the pri- 
mary critical curve of order m = 1 is the classical Ein- 
stein ring. If a source at distance r s lies on the primary 
caustic tube and is large enough that a complete cross- 
section of the caustic tube lies inside the source, then the 
observer will see the primary image and the secondary 
image merged together to cover the whole Einstein ring. 
Since higher order caustics have a greater extension, the 
case that a source intercepts a complete cross-section of 
the caustic tube can only occur for sufficiently low spins 
or very large sources. Therefore, it is much more dif- 
ficult to have a complete Einstein ring of higher order. 
Nevertheless, critical curves are the loci where additional 
pairs of images appear or disappear when a caustic cross- 
ing occurs. Since such images are maximally magnified 
during the process of creation or destruction, it is impor- 
tant to study critical curves because they indicate the 
places where higher order images are most likely to be 
detected. However, the average radius of a critical curve 
depends on the source distance. This is particularly ev- 
ident for the primary critical curve, whose average ra- 
dius at large source distances scales as the Einstein angle 
9e = y/AMr s /r (r + r s ). Higher order critical curves 
tend to an asymptotic shape for large source distances, 
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FIG. 20: Critical curves for a Kerr black hole with spin a = 0.9998A/ and an equatorial observer. The coordinates (0i, #2) in 
the observer's sky are in units of 2M/r a . (a) Primary critical curve for various source distances: starting from the outside, 
r s — 5M, 3M, 1.8A/, 1.1M. (b) Second order critical curve for various source distances: starting from the outside, r s = 
00, 4.2M, 1.6M, 1.1M. (c) Critical curves of order m = 2, 3, 4 for a source at infinity, (d) Critical curves of order m = 1, 2, 3 for 
a source at r 3 = 1.1M. The shadow border is represented by the dashed curve. 



whereas they become smaller for smaller source distances. 
Typically, higher order critical curves are slightly larger 
than the shadow of the black hole, but when the source is 
closer than the unstable photon orbits, the critical curves 
appear entirely inside the shadow of the black hole. 

In Fig. we show some critical curves for a nearly 



extremal Kerr black hole and an equatorial observer. As 
well-known [2|], the shadow is shifted to the right (retro- 
grade side) with respect to the black hole position, with 
the left (prograde) side more flattened. 

In Fig. 120a we show the primary critical curve for var- 
ious source distances. If the source is far from the black 
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hole, the critical curve is also far from the shadow. If 
the source is closer to the black hole, the critical curve 
shrinks gradually and finally enters the shadow when the 
source is closer to the black hole. Note that the critical 
curve enters the shadow from the retrograde side first, 
because the retrograde unstable orbit lies at a larger ra- 
dius. The prograde unstable orbit is closer to the horizon 
and therefore the critical curve stays outside the shadow 
on the prograde side until the source becomes very close 
to the horizon. We can also note that the primary criti- 
cal curve maintains its nearly circular shape at all source 
distances. 

Higher order critical curves, instead, follow the shape 
of the shadow more closely. In Fig. 120b we show the 
second order critical curve for various source distances. 
At very large source distances each critical curve tends 
to a fixed asymptotic curve at finite distance from the 
shadow. As the source distance is decreased, the critical 
curve shrinks and finally enters the shadow. Note that 
higher order critical curves have a small concavity on the 
prograde side, which disappears when the source is closer 
to the black hole. 

In Fig. l2T)b we show the critical curves of order 
to = 2,3,4 altogether for a source at infinity. All criti- 
cal curves are outside the shadow with higher order ones 
closer to the shadow. In Fig. |2"0"B we show the critical 
curves of order to = 1,2,3 for a source at r s = 1.1M 
i.e. very close to the horizon. In this case the primary 
curve is the most internal one, with the higher order ones 
gradually closer to the shadow. 

Finally, the dependence of the critical curves on the 
spin inclination closely follows the dependence of the 
shadow, which becomes more and more circular and sym- 
metric as the observer's latitude increases. 



VII. CONCLUSIONS 

In all applications of gravitational lensing, the study 
of critical curves and caustics of specific lens models has 
always represented a fundamental step in the compre- 
hension of the whole phenomenology. As it can be eas- 
ily imagined, the derivation of caustic surfaces in a full 
general relativistic context is much more involved than 
in classical lens models analyzed under the weak deflec- 
tion paradigm. The simplest general relativistic lens is 
the Schwarzschild black hole. In this model, however, 
the caustics are degenerate and are therefore trivially 
tractable. The first general relativistic lens with a non- 
trivial caustic structure is the Kerr black hole. Going 
beyond the results of previous works, focused on partic- 
ular limits of the caustics 

[H HE EE EE HE] > this paper 
contains a complete investigation of the full caustic struc- 
ture of the Kerr metric. This represents a considerable 
step forward for gravitational lensing phenomenology in 
full General Relativity. 

In summary, we have shown 3-dimensional pictures of 
the primary and higher order caustic surfaces. We have 



analyzed the dependence of the cross-sections of the caus- 
tic surfaces on the source distance, the black hole spin, 
the spin inclination and the caustic order. We have shown 
that the caustic surfaces always wind an infinite number 
of times around the horizon following a logarithmic law. 
The size of the caustics always remains finite in the whole 
parameter space and there are no transitions to different 
kinds of caustic singularities in the lens mapping. For 
extremal spin values, the size of higher order caustics in- 
creases exponentially with the caustic order. This implies 
that the number of higher order images also grows expo- 
nentially, as the same caustic can be crossed an exponen- 
tial number of times. We have compared our results with 
previous analytical and numerical approximations, find- 
ing perfect agreement in the respective limits of validity 
of the approximations. We have also shown some critical 
curves in the observer's sky, focusing on their dependence 
on the source distance. 

In addition, the code developed in this paper for the 
calculation of the caustics is particularly well-suited for 
the study of higher order images, as being inspired by 
the Strong Deflection Limit methodology. An interesting 
future development of this code might be the implemen- 
tation of an efficient resolution algorithm of the Kerr lens 
equation including higher order images. 

Besides the purely theoretical interest of the analysis 
presented in this paper, we can easily imagine that the 
knowledge of the complete caustic structure of the Kerr 
black hole in all ranges of parameters will stand as an ex- 
tremely helpful guide in future astrophysical applications 
involving very strong light bending by the gravitational 
field of spinning black holes. 
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APPENDIX A: FROM THE GEODESICS 
EQUATION TO THE LENS MAPPING 

In this appendix, we describe all the steps of the cal- 
culation of the caustics in Kerr spacetime. Section IA II 
recalls the basic facts about Kerr geodesies. Section [A 21 
discusses the unstable photon orbits and the choice of 
good variables for numerical integration. Section IA 31 
deals with angular integrals, Section [A 41 with radial in- 
tegrals. Section [A 51 introduces the ip variable and shows 
the lens mapping explicitly. Finally, Section fA 6l discusses 
details on the numerical implementation of the strategy 
for searching critical points and drawing caustic surfaces. 
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1. Basics of Kerr geodesies 

Thanks to the separability of the Hamilton- Jacobi 
equation [|3] , the geodesic motion for a massless particle 
is described by two integral equations 



± 



dr 



= ± 



aJ 



csc 2 $ 



dx — a 



d.r 



(Al) 
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where 

e 
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a 2 cos 2 ■& - J 2 cot 2 i? 



(A3) 
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J and Q are two constants of motion. The first can be 
identified with the projection of the angular momentum 
of the particle on the spin axis and the second is related to 
the square of the total angular momentum. The double 
signs before the integrals are chosen in such a way that all 
pieces of the integrals between two consecutive inversion 
points are summed with a positive sign. 

A photon reaching the observer with constants of mo- 
tion J and Q will then be detected by the observer at 
angular coordinates 



J 



r oA /l-Mo' 



' + K K - 



1-M 2 
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Only an indetermination on the sign of 82 remains. 

Conversely, the constants of motion of a photon reach- 
ing the observer from angles (8i, 82) are found by invert- 
ing Eqs. jM])-CM} 



J=-9 iro ^T 
Q = 6 2 r 2 +£(6 2 r 2 -a 2 ). 
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2. Choosing good labels for photon trajectories 

All photon trajectories ending at the observer can be 
completely specified by assigning proper values to the 
constants of motion J and Q. Equivalently, one can 
specify the coordinates in the observer sky 8\ and 82- 
However, in order to keep numerical errors low, we may 
need to find new labels for the photon trajectories. In the 
case of the Kerr metric, the most challenging geodesies 
are those approaching the unstable photon orbit. We 
have found that using variables inspired by the Strong 



Deflection Limit (SDL) technique 19] proves very conve- 
nient, opening the way to a reliable description of higher 
order caustics. Therefore, in this subsection we recall the 
definition of the SDL parameters S and e, and introduce 
the variable 77 which replaces the £ of Ref. [19( . 

The inversion points in the radial motion are obtained 
by solving Eq. R = 0. Photons reaching an observer at 
infinity can have one or zero radial inversion points. If 
dR/dr — is satisfied at the same time with R = for 
some value r = r m , then the photon remains at radial co- 
ordinate r m forever. However, small perturbations would 
make the photon fall into the black hole or escape to in- 
finity, because the circular photon orbits around Kerr 
black holes are unstable. The unstable photon orbits 
are characterized by particular values of the constants of 
motion J and Q. In order to find them, we can solve 
equations R — and dR/dr — for J and Q as func- 
tions of r = r m . As r m varies, J and Q describe a locus 
in the (J, Q) space, corresponding to all possible unstable 
photon orbits. Explicitly, this locus is given by Q 



^m(^m) 
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A photon orbiting at r m is thus characterized by con- 
stants J m (r m ) and Q m (r m ). The allowed values of r m 
are those keeping Q positive. By perturbing the unsta- 
ble orbit, we can make the photon escape and reach a 
distant observer. In this case, plugging Eqs. (|A9[) - (|A10[) 
into Eqs. (|A5|) - (|A6|) . one finds that the escaped photon 
is detected at position 



^l.m (^*m ) 

h,m{r m ) 
A(r m ) 
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(A12) 

2 

°(A13) 



Not all perturbed photons can reach an observer out of 
the equatorial plane. Only those photons with r m such 
that A(r m ) > can be detected. Solving Eq. A(r m ) = 
numerically, we can find the two extrema of the allowed 
interval for r m for any given values of the parameters a 
and [i - We shall indicate these two extrema by r+ and 
r_. 

Letting r m vary between r + and r_ in Eqs. (jAlip and 
(|A12|) and allowing both signs for 82, we obtain a closed 
curve in the observer's sky, which is usually referred to 
as the border of the "shadow" of the black hole, because 
the radiation deflected from the black hole would appear 
outside this closed curve. Some shadows borders for dif- 
ferent valus of a are shown in Fig. [TJ When the observer 
lies on the equatorial plane (jj, = 0), r + coincides with 
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the radius of the unstable photon orbit for equatorial 
prograde photons, whereas r_ becomes the radius of the 
unstable photon orbit for equatorial retrograde photons, 
which is typically larger. Increasing the inclination, r + 
and r_ approach each other until they coincide when the 
observer lies at the pole of the black hole (/i = 1). 

As r+ and r_ depend on the spin and the inclination, 
it is convenient to introduce the variable 77, related to r m 
by 



become 



1 



r_l_(l — cos 77) + r_ (1 + cos 77)] 



(A14) 



As 77 varies from to tt, r m varies from r_ to r + . If we 
replace the double sign in Eq. (|A12|) by sign [77], we can 
get the whole shadow border at once, by varying 77 in the 
range [— tt, tt]. 

As anticipated in Section II, 77 works as an angular 
coordinate: 77 = corresponds to retrograde photons, 
appearing on the right of the black hole; 77 — ir/2 cor- 
responds to photons on nearly polar orbits, reaching the 
observer from above the black hole; 77 = tt corresponds to 
prograde photons, detected on the left of the black hole. 

In addition, we introduce the new variable e by the 
relations 



h =0i, m (7?)(l + e) 
h = e 2 , m {ri){l + e). 



(A15) 
(A16) 



As e varies in the range [—1, +00) and 77 in [— tt, tt], B\ and 
62 span the whole observer sky. This construction is very 
close to that presented in Ref . [l9[ , though here we make 
no expansion in powers of a. The variables 77 and e can 
be used to label geodesies ending at the observer and are 
particularly well-suited for numerical calculations. The 
corresponding values of the constants of motion J and 
Q can be found by plugging Eqs. (|A15[) and (|A16|) into 
Eqs. CEO) and 

When e > 0, Eq. R = admits a non-degenerate 
inversion point 7%, whose expression can be written as 



r* = r m (l + 5). 



(A17) 



S is a function of e and 77. It can be found by solving 
Eq. R = numerically for given values of e, 77 and the 
parameters a and [i Q . For small values of e, we have 
6 ~ \ft. However, for low values of the inclination fi a and 
for high values of the spin, prograde photons (77 ~ tt) tend 
to satisfy a linear relation of the type S ~ e. However, 
S is not an observable and these relations are specific of 
Boyer-Lindquist coordinates. 



3. Angular integrals 

Eqs. (|A1[) and (|A2|) contain integrals on the polar an- 
gle 7?. Performing the transformation cos § = fi, they 



Ji = ± 




where 
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The integration extrema are the source and observer 
polar coordinates fj, s and fi - However, for each inversion 
in the polar motion, we must split the integration domain 
and change sign. 

The final result can be expressed in terms of elliptic 
integrals of the first and third kind 
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where m is the number of inversions in the polar motion, 
and the sign[?7] takes into account the fact that for posi- 
tive 77 the observer is reached from above and for negative 
77 is reached from below. We recall the definitions of the 
elliptic integrals 



F(X,k) 



(It) 
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U(n,k) = U(n,w/2,k). 



(A31) 



Besides depending on the source polar coordinate /i s , 
the functions J\ and Ji also depend on e (weakly) and 77 
through the constants of motion J and Q, appearing in 
the expressions of fi+ and /z_ . 
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4. Radial integrals 

The geodesies equations (|A1|) and (|A2[) also contain 
two integrals on the radial coordinate 



h = 



h = 



dx 

Vr 

dx 
1 



a J 



aVr 
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(A32) 
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These integrals depend on the labels (e, rf) identifying 
the geodesic through J and Q. 

The photon might be emitted from the source either 
outward or inward. We must therefore distinguish these 
two cases and solve the integrals accordingly. 

If the photon is emitted in the outward direction, there 
is no inversion point and the integrations are carried out 
in the domain [r s ,oo). For numerical reasons, it is con- 
venient to change the integration variable to z, defined 
by 

r = (A34) 
1 — z 

Then the integration range becomes [1 — r m /r Sl l]. In 
this case, the integrand is sharply peaked at z — [201 ]. 

If the photon is emitted inward, it must have an inver- 
sion point r* somewhere between r m and r s in order to 
come back and reach a distant observer. Therefore, this 
case is relevant for lensing only if e > and r s > r m . 
The radial integrals are split in two pieces: the approach 
piece, with extrema [r*,r s ], and the escape piece, with 
cxtrcma [r».oo). The two pieces arc conveniently calcu- 
lated using the variable 

(A35) 



l-z' 

which maps the integration domains to the finite inter- 
vals [0, 1 — r*/r s ] and [0, 1] respectively. In this case the 
integrand diverges at z = 0, but the result of the inte- 
gration is finite for any value of e > 0. As e — > 0, the 
integrals diverge logarithmically, except for prograde or- 
bits at high values of the spin, for which the integrals 
diverge as e . 

For a given source at r s > r m , there exists a maximum 
value for e (we call it e max ) such that R is positive only 
for e < e max - When e = e max , r s becomes a solution of 
the equation R = and thus coincides with the inversion 
point r*. 

Summing up, the functions Ii{e,rj) and ^(e,??) have 
two branches if r s > r m : one corresponding to photons 
emitted outward, which exists for e in [—1, e maa ;], and the 
other corresponding to photons emitted inward, existing 
for e in (0, e max ] . The two branches nicely join at e max . If 
r s < r m , only the outward branch exists for e in [—1,0). 



5. Lens mapping 

At this point, we introduce a new variable tp, which 
considerably simplifies numerical calculations and as- 




FIG. 21: Relation between tp and e for a = 0.2M, {i = 0, 
77 = 7r/4. From left to right, the curves are drawn for r s = 
2M, 3M, 4M, 5M, 6M respectively. Note that for these values 
we have r m — 3.16M. 



sumes a key role in their success. This variable is defined 
by the relation 



2K{k)i) = afx-h(e, n) + sign[77]F(A , k). 



(A36) 



The variable ip is directly defined in terms of the ra- 
dial integral 7i, thus eliminating the problems of the 
two branches that this function has when is expressed 
in terms of e. The reason behind the coefficient 2K(k) 
will be clear later. 

In Fig. [21] we plot ip vs e for a generic choice of the 
parameters. We note that when r s < r m , e < for all 
values of 'ijj, whereas for r s > r m e reaches the maximum 
value e max , corresponding to r s — and then tends to 
zero for very large values of ip. On the other hand, e is a 
single- valued function of tp. Therefore, we can replace all 
occurrences of e by e(ip) obtained by solving Eq. (|A36j) 
numerically. 

Multiplying Eq. (|Al|l by a/i_ , separating the /z s term 
and applying the Jacobi elliptic function sn, we get 



[i s = sign[?7]^ + sn [2K(k)ip, k] 



(A37) 



which represents the lens equation for the polar coordi- 
nate fi s . 

The Jacobi elliptic function sn(a;, k) is periodic in x 
with period 4K(k). It tends to the usual trigonomet- 
ric sin when k is small (which is more or less our case, 
since fi + < fx— ). So, the right hand side of Eq. (|A37|) 
vanishes at if) = m with m integer, is maximum at 
tp = 2m + 1/2, and minimum at tp = 2m + 3/2. The 
interval in which there are m inversion points in the po- 
lar motion is m — 1/2 < tp < m + 1/2. As we shall see, 
the fact that this interval has a trivial form in terms of tp 
for any values of m is crucial in the search for the critical 
points of higher order. If we had kept e, the definition of 
the search domain would have been much more problem- 
atic. In Ref. [l^], the variable tp played a similar role, 
but its definition was restricted to the limit of small black 
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hole spin. In that case ip reduces to the equivalent de- 
flection by a Schwarzschild black hole of the same mass. 
Eq. (|A36j) generalizes the former definition and allows 
a fruitful use of this variable in all the Kerr parameter 
space. 

As regards the azimuthal motion, the extrema of <f> are 
4>i — 4> s and <pf = <j) . Recalling that we have set <p = 0, 
Eq. (|A2[) can be written in the form 

S =a[h(i>,T 1 ) - J 2 (iM)] ~ JJ 2 (a*.,iM). (A38) 

This is the lens equation for the azimuthal coordinate <p s . 
Along with (1A37|) . it constitutes the lens mapping for a 
Kerr black hole in the form of Eqs. © and 0. 

6. Remarks on numerical implementation 

Here we give a few technical details regarding the nu- 
merical implementation of the algorithm for finding crit- 
ical curves and caustics. 

As input parameters we consider the black hole spin 
a, the inclination fi , the source distance r s , the caus- 
tic order m and the value of the variable rj. Of course, 
depending on which kind of plot we want to draw, we 
may decide to cycle on different variables. For example, 
by cycling on ?y keeping all other parameters fixed, we 
get a cross-section of the caustic surface at fixed source 
distance. Cycling on a, while keeping rj at some specific 
value, we may get the position of a cusp (for example) as 
a function of the black hole spin. 

Starting from each set of parameters {a, ji a , r s , to, 77}, 
first we calculate r + and r_ by solving Eq. A(r) = 



with A(r) given by Eq. (|A13|) . Then we find r m by Eq. 
(IAT41) . 

For any test value of we numerically invert Eq. 
()A36[) to find e, keeping track of the branch of Ji in 
which we are. If we are in the branch characterized by 
the presence of the radial inversion point, we also need 
to calculate 5 by numerical solution of Eq. R = 0. 

All these equations are solved using the Mathematica 
FindRoot routine with the secant method. The integra- 
tions of I\ and I2 are performed using the NIntegrate 
routine, which works pretty well on the finite domains 
obtained by changing to the z integration variable. 

At this point, we can evaluate the lens equations (|A37|) 
and (|A38[) to get fi s and cf> s , without further effort. 

The Jacobian of the lens mapping is obtained by eval- 
uating difference quotients of ji s and <j) s with respect to 
ip and 77. 

Finally, the zero of the Jacobian at fixed rj is found 
by applying the secant method with respect to ip m t he- 
range [to — 1/2, to + 1/2] . 

Depending on the required precision and the number 
of points in the cycle, the evaluation of a caustic may 
take several minutes. We have required a precision of 
10~ 6 in tj] and about 500 points for each caustic, refining 
the sampling in the neighborhood of the cusps. With this 
setup, a cross-section at fixed r s takes about 5 minutes on 
a laptop. Numerical noise has been detected on prograde 
orbits at extremal spin only starting from the fifth order 
caustic. 
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